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Abstract 

The  effect  of  incoming  organized  disturbance  and  free-stream  turbulence  on  leading-edge  heat  transfer  is 
investigated  numerically.  An  optimum  length  scale  is  found  to  give  the  maximum  heat  transfer  enhancement 
for  the  organized  disturbance  case.  Beyond  this  optimum  value,  the  enhancement  decreases  with  the  increase 
of  length  scale.  For  the  free-steam  turbulence  case,  large  eddy  simulation  with  dynamic  SGS  model  is 
performed  at  Reynolds  number  Reo  =  104  based  on  upstream  velocity  u0 0  and  the  leading  edge  diameter 
of  curvature  D.  The  free-stream  turbulence  is  specified  as  homogeneous,  isotropic  turbulence  with  intensity 
u'ma/u oo  =  0.08  and  integral  length  scale  L/D  =  0.1.  Three  different  regions  characterize  the  interaction 
of  turbulence  impinging  on  the  leading  edge.  For  the  conditions  of  the  simulations  a  turbulent  heat  transfer 
enhancement  of  11%  is  obtained,  which  is  in  fair  agreement  with  the  experimental  data.  The  level  of  heat 
transfer  enhancement  is  modest  because  of  the  Reynolds  number  is  low.  However,  our  results  show  that  in 
the  region  very  close  to  the  leading  edge,  the  energetic  turbulence  length  scale  decreases  to  the  order  of  2-3 
times  the  local  boundary  layer  thickness,  so  a  high  grid  resolution  is  needed  for  accurate  prediction  of  heat 
transfer  using  large  eddy  simulation.  This  is  a  challenge  for  future  investigations  where  simulations  at  higher 
Reynolds  numbers  and  transonic  flow  conditions  are  planned. 

Our  results  motivate  a  hybrid  simulation  strategy  where  the  tui  oulence  outside  and  away  from  a  blade 
surface  is  captured  using  LES  techniques  while  a  finer  DNS-like  grid  is  embedded  within  the  near-wall  region 
to  resolve  the  smaller  eddies  responsible  for  near-wall  effects.  Such  a  methodology  is  being  developed  in  an 
extension  of  the  work  supported  under  this  grant. 

1  Introduction 

The  effect  pf  significant  free-stream  turbulence  generated  by  the  combustor  and  the  upstream  blade  rows 
on  turbine  blade  heat  transfer  is  a  critical  problem  in  advanced  gas  turbine  design.  An  improved  under¬ 
standing  and  accurate  prediction  of  the  heat  transfer  under  free-stream  turbulence  is  of  essential  importance 
for  achieving  high  turbine  inlet  temperatures  and  better  thermal  efficiency.  It  is  known  that  free-stream 
turbulence  can  cause  a  large  increase  in  leading  edge  heat  transfer.  Many  experiments  have  been  conducted 
to  investigate  the  correlation  of  the  heat  transfer  enhancement  to  various  flow  parameters,  such  as  Reynolds 
number,  Mach  number,  turbulence  intensity,  turbulence  length  scale  and  the  leading  edge  geometry  [1-5]. 
Ames  et  al.  [4]  found  that  both  turbulence  intensity  and  turbulence  length  scale  are  important  in  correlating 
stagnation  region  heat  transfer.  Van  Fossen  et  al.  [5]  found  the  heat  transfer  augmentation  increases  with 
decreasing  length  scale  and  the  spanwise  averaged  heat  transfer  coefficient  can  be  more  than  50  percent 
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higher  than  the  two-dimensional  laminar  value.  To  understand  the  sensitivity  of  the  heat  transfer  in  the 
stagnation-region  boundary  layer  to  free-stream  disturbance,  Sutera[6]  proposed  a  theoretical  model  to  study 
the  vorticity  amplification  due  to  the  stretching  of  the  vortices  in  the  diverging  stagnation-flow.  Rigby  et 
al  [7]  numerically  simulated  spanwise  variations  in  the  laminar  free-stream  in  the  flow  over  circular  and 
elliptical  leading  edges  and  found  the  wall  heat  transfer  rate  increased  by  as  much  as  25  percent.  Bae  et  al [8] 
performed  a  numerical  simulation  of  laminar  plane  stagnation-region  flow  with  organized  free-stream  span- 
wise  disturbances,  and  found  three  different  regimes  for  the  amplification/damping  of  streamwise  vortices 
depending  on  the  disturbance  length  scale. 

Numerical  simulations  using  organized  laminar  disturbances  can  isolate  various  factors  such  as  flow 
disturbance  intensity  and  length  scale  that  affect  leading  edge  heat  transfer,  and  provide  a  simple  but 
approximate  analogy  to  the  unsteady  interacting  process  between  free-stream  turbulence  anu  the  leading 
edge.  However,  the  direct  computation  of  free-stream  turbulence  interacting  with  the  leading  edge  is  very 
important  for  understanding  the  global  mechanism  of  turbulent  heat  transfer  enhancement,  and  for  further 
developing  physically-based  turbulence  models  capable  of  accounting  for  the  effects  of  free-stream  turbulence. 
In  the  present  study,  we  perform  a  large  eddy  simulation  (LES)  of  free-stream  turbulence  interacting  with 
the  elliptical  leading  edge  of  a  finite  thickness  isothermal  flat  plate.  The  flow  configuration  is  based  on  the 
experiments  previously  studied  by  Van  Fossenef  al.[ 5].  The  purpose  is  to  explore  in  detail  the  unsteady 
interaction  process  between  free-stream  turbulence  and  leading-edge  flow^  and  obtain  a  better  understanding 
of  the  heat  transfer  augmentation  mechanism.  In  addition,  the  simulation  also  serves  as  an  evaluation  of  the 
feasibility  of  using  LES  as  a  tool  for  accurate  heat  transfer  prediction  in  the  gas  turbine  environment. 

The  report  is  organized  as  follows.  The  governing  equations  and  methodology  for  large  eddy  simulation 
are  described  first,  followed  by  the  numerical  schemes  for  the  implicit  time  marching  and  high  order  spatial 
differencing.  The  generation  of  free  stream  turbulence  with  specified  intensity  and  length  scale  is  discussed 
together  with  the  inflow^  boundary  conditions  for  specifying  the  free-stream  turbulence.  Both  the  results  from 
laminar  organized  disturbances  and  free-stream  turbulence  are  presented  and  discussed.  The  heat  transfer 
enhancement  results  are  compared  with  reported  experimental  values,  and  the  typical  flow  structures  at  the 
leading  edge  are  investigated.  Concluding  remarks  summarize  the  results  and  discuss  future  plans  for  this 
research. 


2  Governing  Equations 


The  governing  equations  for  large  eddy  simulation  are  the  Favre-filtered  compressible  continuity  and  Navier- 
Stokes  equations,  written  in  the  non-dimensional  primitive  variables: 


P,t  +  {pui),i  =  0 


PUi,t  +  PUjUiJ 
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with  the  filtered  equation  of  state  for  ideal  gas 


(3) 


(4) 


In  the  above  equations,  p  is  the  density,  u2  is  the  velocity  vector,  p  is  thermodynamic  pressure.  M,  Re  and 
Pr  are  the  Mach  number,  Reynolds  number,  and  Prandtl  number  defined  as 


M  _  Woo  Re  _  PocUccD_  Pr  =  Moo Cp 

Coo  Moo  K'OO 


(5) 


Subscript  oo  denotes  the  far  upstream  ambient  conditions,  c  is  the  speed  of  sound.  D  is  the  leading 
edge  diameter  of  curvature,  p  and  A  are  the  first  and  second  coefficients  of  viscosity,  and  k  is  the  thermal 
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conductivity.  50  =  +  ujfi)  is  the  rate-of-strain  tensor.  Subscripts  following  a  comma  denote  partial 

differentiation  with  respect  to  the  subscript,  and  the  Einstein  summation  convention  is  used.  The  filtered  or 
large  scale  flow  quantity  /,  denoted  by  an  overbar  is  given  by: 

7=  [  G(x-x')f(x’)dx'  (6) 

Jd 

where  G  is  some  spatial  filter-kernel  and  the  integral  is  over  the  flow  domain.  The  Favre-averaged  quantity, 
denoted  by  a  tilde  is  defined  as 

f  =  pf/p  (?) 

The  effect  of  the  small  scales  are  present  in  the  above  equation  through  the  sub-grid  scale  (SGS)  stress  tensor 
Tij  in  the  momentum  equation, 

Tij  =  p(UiUj  -  UiUj)  (8) 

and  the  SGS  heat  flux  in  the  energy  equation, 

qt  =  p(uiT  -  inf)  (9) 

Several  other  subgrid-scale  terms  have  been  neglected  from  the  energy  equation  following  the  same  argument 
as  in  [9].  The  power  law  for  the  molecular  viscosity  coefficient  p  is  employed 

p  =  f07  (10) 

and  Stokes’  hypothesis  is  used  for  the  second  viscosity  coefficient  A. 


3  Dynamic  SGS  model 

To  close  the  momentum  and  energy  equation,  Tij  and  qi  must  be  modeled.  For  the  trace-free  Smagorinsky 
eddy  viscosity  model  is  used  - 

Tij  =  ±q%  -  2CpA2|5|(Sy  -  ±S**<5tj)  (11) 

where  q 2  =  t„  is  the  isotropic  part  of  the  SGS  Reynolds  stress  tensor,  Sij  =  ^{uij  4-  Uj,i)  and  \S\  = 
(2 SijSij)1/2.  The  SGS  energy  q2  is  parameterized  using  Yoshizawa’s  expression: 

■  q2  =  2C7/?A2|S|2  (12) 


For  the  qi ,  the  eddy  diffusivity  model  is  used 


Qi  =  - 


CpA2\S\ - 
1 


(13) 


where  Prt  is  the  SGS  turbulent  Prandtl  number  and  C  is  the  same  as  in  (11).  We  use  the  dynamic  procedure 
to  compute  the  eddy  coefficients  C,  C7  and  SGS  turbulent  Prandtl  number  Prt ,  which  are  all  functions  of 
instantaneous  flow  variables.  The  key  concept  of  the  dynamic  model  is  the  effective  utilization  of  the  large- 
scale  field  which  is  computed  directly.  This  information  is  brought  to  bear  by  introducing  a  “test  filter” 
with  a  larger  filter  width  than  the  resolved  grid  filter,  which  generates  a  second  field  with  scales  larger  than 
than  the  resolved  field.  Details  of  the  derivation  may  be  found  in  [9].  Let  a  spatially  test-filtered  quantity 
be  denoted  by  a  caret  and  the  width  of  the  test  filter  is  denoted  by  A,  the  SGS  coefficients  (7/,  C  can  be 
computed  as 


where 


c,=  <  Ckk 

(14) 

<  2(£A2|S|2  -  A2p|S|2)  > 

n  _  <  (£*j  —  3  £kk$ij)Mij  > 
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and 


■Mij  —  —  2pA“|S|(5ij  —  ^Skk^ij)  + 

2A2(p\5\'sij  -  ±p\§\Skk6ij)  (17) 

The  turbulent  Prandtl  number  Prt  is  given  by 

Prt  =  c  <  -  ggjji Ma  >  (18) 

<  (  rhout  pT  /p-  puiT)f  i  > 

where  C  is  given  by  (15),  and  <>  indicates  a  volume  averaging  procedure  which  is  needed  to  make  the 
determination  of  C\  and  the  other  SGS  coefficients  well  conditioned.  For  all  of  the  computations  presented 
here,  the  ratio  of  the  test-filter  width  to  the  grid-filter  width  A/ A  is  taken  as  2.0  and  Cj  =  0.0. 


4  Numerical  Method 

The  solutions  to  the  Favre-averaged  equations  (l)-(4)  are  obtained  numerically  using  implicit  time  integration 
with  an  approximately-factorized  finite  difference  scheme.  In  the  case  of  time  accurate  computations,  the 
Newton-like  sub-iteration  scheme  is  employed.  For  the  purpose  of  analysis,  the  Navier-Stokes  equations  are 
first  cast  into  a  general  form: 

U,t  +  F(U)  =  0  (19) 

where  U  is  the  vector  of  flow  variables  and  F(U)  represents  the  nonlinear  and  viscous  terms.  An  implicit 
approximation  in  time  for  this  general  equation  can  be  written  as  [10] 

[/  +  ^^-A(Un)]AUn  =  --^—F(Un)  +  — ^ — At/”-1  (20) 

1  1  +  v?  1  -t-  1  +  <f 

where  A  =  F(U)tu  is  the  Jacobian  of  F.  This  equation  is  kept  in  what  is  called  “Delta  Form”,  referring  to 
the  A Un  =  Un+1  —  Un  term  on  the  left  hand  side  and  is  second  order  accurate  in  time  for  #  =  =  0  or 

It  forms  the  basis  for  most  implicit  time  integration  schemes. 

To  represent  the  complex  geometry  using  a  body-fitted  mesh  and  clustering  points  in  the  regions  of  high 
spatial  gradient  such  as  the  boundary  layer,  a  global,  two-dimensional  mapping  is  constructed  from  physical 
space  (x,y,z)  to  a  uniform,  Cartesian  computational  space  (£,77,  z).  Both  £  and  rj  are  in  the  range  [0,1].  In 
computational  space  (f,?7,z),the  Jacobian  A(U)  can  be  written  as  [11] 


A(U)  =  AU't  +  BUr,  +  CUxZ  +  DU 
-VuUM-VtaUtn-VisVts 

-VnflUm  -  Vr)zU,T)z  ~  VzzUlZz  (21) 

where  U  =  {p,u,v,w,T}T  is  an  unsteady  flow  field  and  the  matrices  (A,  F,C,  D,Vij)  are  function  of  U  and 
its  gradients  and  also  contain  the  metrics  terms. 

Substituting  (21)  into  (20),  and  using  conventional  finite  difference  or  finite  volume  schemes  for  the  spatial 
derivatives  usually  results  in  large  sparse  matrix  systems,  which  are  not  too  difficult,  but  prohibitively  time 
consuming,  to  solve  directly.  Approximate  factorizations  are  employed  to  improve  the  efficiency,  reduce  the 
computer  storage  requirements  and  enhance  numerical  stability.  Thus,  factorized  in  the  f,  r\,z  directions, 
equation  (20)  can  be  written  as: 


[/  +  aAt{A&t  +  D-  VUA a)]P  =  K 

(22) 

[I  +  aAt(BAn-VnriAnn)}Q  =  V 

(23) 

[I  +  aAt(CAz  -  V..A«)]A Un  =  Q 

(24) 
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where 


(25) 


a  = 
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n  =  - 
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1  +  f 


F(Un)  + 
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A  U 
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The  factorizations  introduce  a  0(At~)  error  compared  to  the  unfactorized  equations.  This  factorization  error 
can  be  eliminated  by  employing  a  subiteration  time  stepping  procedure  to  recover  the  solution  of  (20).  For 
this  purpose,  rewrite  equation  (20)  by  substituting  a  new  iterative  index  p+  1  for  the  n  +  1  terms  and  adding 
Up  to  both  sides  of  the  equation,  to  obtain 


[/  +  A(UP)]AUP  =  —£Lf(Up)+ 

1  +  if  l+ip 

( Up  -  +^-Un  +  —¥—Un~l) 

v  1  +  <p  l  +  ip 


(26) 


with  AUP  =  ( Up+l  -  Up )  and  A{UP)  =  dF/dUp.  Equation  (26)  is  the  basic  sub-iteration  time  advancement 
scheme  which  will  yield  second  order  time  accuracy  independent  of  choice  of  A(UP)  if  the  sub-iteration 
process  converges.  The  n  and  n  —  1  terms  are  evaluated  from  previous  time  level  and  after  iterating  p  times 
the  solution  at  time  level  n- hi  will  be  taken  from  the  most  recent  Up+l.  For  example,  choosing  p  =  |  and 
d  =  1,  in  the  limit  (assuming  the  iterative  process  converges)  UpJrl  =  Up ,  setting  t/n+1  =  Up+1 ,  we  have 


3 £/n+1  ~4 Un  +Un-X 

2At 


+  F(Un+  )  =0 


(27) 


which  is  a  second  order  in  time,  fully  implicit  approximation  of  equation  (19). 

This  sub-iteration  scheme  requires  that  first  the  iterative  process  converge,  and  second  the  local  sub- 
iteration  process  be  performed  until  the  error  is  of  higher  order  in  time.  In  practice  one  finds  that  two  to 
three  sub-iteration  gives  adequate  accuracy  for  most  problems. 

For  the  spatial  discretization,  consider  derivative  operators  in  the  uniform  computational  space  (£,r/) 
where  £i,r)i  represent  the  nodal  locations  with  &  =  A£(z  -  1)  for  1  <  i  <  and  r]j  ~  Ar](j  -  1)  for 
1  <  j  <  At  interior  nodes,  the  fourth-order-accurate  central  difference  scheme  is  used  for  first  and 
second  derivatives. 

i2^[(/i— 2)  -  ft+ 2)  -  8(/i-i  -  fi+ 1)]  (28) 


d-l) 


(|0  — L 


12(A#[_(/i_2"2/,+/<+2 

+  16(/i-i  -  2 fi  +  /t+i)] 


(29) 


Near  the  computational  boundaries,  finite  difference  operators  that  are  biased  toward  the  interior  are  re¬ 
quired.  As  for  the  interior  points,  a  five  point  stencil  is  used  at  the  boundary  grid  point  .  The  resulting 
difference  schemes  are  forth  and  third  order  accurate  for  the  first  and  second  derivatives  respectively.  Hence, 
At  the  first  grid  point  £i ,  the  first  and  second  derivative  are  expressed  as  : 


{%) ,  =  l2^[_25/l  +  48/2 ' 36/3  +  16/4  ~  3/s! 
{w) ,  =  J2^[Ufl  ~  20/2  +  6/3  +  4/4  "  h] 

while  at  the  second  point  £2,  they  are 

'|)3  =  TSf[-3/,-10/2  +  ‘Ws-6/<+/sl 

77[35/i  -  104/a  +  ll/s  -  56/4  +  ll/s] 


U ? 


1 

12(A0'2 


(30) 

(31) 

(32) 

(33) 
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Figure  1:  Flow  configuration  schematic 


Similar  expressions  hold  for  the  derivatives  at  nodes  of  -  1  and  but  with  the  stencils  reversed  and  the 
signs  switched  on  the  coefficients  for  the  first  derivatives. 

For  the  test  filter  in  the  dynamic  SGS  model,  the  following  Simpson’s  rule  is  used: 

fi  =  g(/i- 1  +  4  fi  +  /*+i)  (34) 

It  is  only  applied  in  the  homogeneous  direction,  where  the  grid  is  uniform  so  that  the  filtering  and  differential 
operation  commute. 


5  Initial  and  Boundary  Conditions 


The  geometry  and  flow  configuration  are  based  on  the  experiments  previously  studied  by  Van  Fossen  et  al.  [5] 
and  is  shown  in  Fig.  1.  A  finite-thickness  fiat  plate  is  placed  horizontally  in  the  center  of  the  test  section 
of  a  wind  tunnel.  The  flat  plate  has  an  elliptical  leading  edge  with  a  major  to  minor  axis  ratio  3:1  and  is 
held  at  a  uniform  temperature.  The  flow  is  confined  on  four  sides  by  walls  which  define  the  test  section. 
The  computational  domain  is  15D,10D  and  ID  in  length,  height  and  width.  The  inflow  boundary  is  placed 
at  4 D  upstream  of  the  leading  edge  and  the  isotropic  homogeneous  free-stream  turbulence  is  convected  into 
the  domain  from  the  inflow  boundary  and  impinges  on  the  leading-edge. 

The  boundary  conditions  used  are  as  follows.  On  the  body  surface,  no-slip  and  isothermal  conditions 
are  applied.  The  density  is  obtained  from  the  continuity  equation.  At  the  outflow,  the  parabolized  Navier- 
Stokes  equation  is  solved  using  a  known  pressure  gradient  taken  from  previously  computed  potential  flow 
solution[ll].  On  the  inflow,  the  boundary  conditions  are  based  on  the  locally  one-dimensional  Riemann 
invariants.  For  a  subsonic  flow,  four  quantities  must  be  specified.  Here  we  constrain  the  entropy,  spanwise  and 
tangential  velocity,  and  the  incoming  Riemann  invariant.  The  locally  one-dimensional  Riemann  invariants 
are  defined  in  the  direction  normal  to  the  inflow  boundary  as 


Ri  =  vn 


2c 

7-1’ 


Ri  —  vn  T 


2c 

7  “  1 


(35) 


where  c  is  the  local  speed  of  sound  and  vn  is  the  velocity  normal  to  the  boundary.  Then  on  the  inflow 
boundary, 


s  =  Soo,  w  =  w‘ ,  vt  =  vt  +  v[ , 

R\  =  R\  +  R\ ,  R‘2  =  R‘2int  (36) 

here  the  over  line  quantities  denote  the  base  flow  in  the  absence  of  inflow  turbulence,  primed  quantities 
represent  the  free-stream  turbulence,  and  the  subscript  t  denotes  the  tangential  velocity  along  the  inflow 
boundary  in  the  x  —  y  plane.  In  the  present  computation,  the  inflow  is  not  far  away  from  the  body  and  the 
flow  passage  is  significantly  blocked  by  the  presence  of  the  flat  plate,  so  the  potential  solution  should  not 
be  used  directly  to  form  the  incoming  Riemann  invariant  due  to  the  development  of  the  boundary  layer  on 
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the  body.  Hence  we  first  compute  the  two  dimensional  steady  base  flow  using  the  same  set  of  the  inflow 
conditions  but  in  the  absence  of  turbulence.  The  velocities  vt  and  vn  are  determined  by 

—  0,  Vn  =  Vnp  (3f) 

Here  ui  is  the  vorticity  at  the  inflow  and  vnp  denote  the  normal  velocity  obtained  from  the  potential  solution. 
These  overline  quantities  are  used  to  form  the  Ri  and  the  primed  quantities  are  interpolated  from  the 
precomputed  isotropic  homogeneous  turbulence.  The  interpolation  is  implemented  using  fourth  order  B- 
spline.  The  /?2m<are  computed  from  first  order  extrapolation  from  the  interior  domain  Him!  —  2/?>,v-i  - 
R-,n_2.  When  the  Riemann  invariants  are  obtained,  the  values  of  vn  and  c  at  the  inflow  can  be  obtained. 
When  passing  the  free  stream  disturbance,  the  turbulence  is  multiplied  by  a  function  exp(-(x/ r>)'2)  so  that 
turbulence  only  enters  the  center  vertical  portion  near  x  =  0  of  the  inflow  plane.  This  is  reasonable  since  the 
turbulence  far  from  the  body  in  the  transverse  direction  is  not  expected  to  impact  the  flow  near  the  body. 

6  Generation  of  Free-Stream  Turbulence 

In  the  present  study,  the  mean  flow  at  the  inlet  boundary  is  nearly  uniform  and  shear-free,  so  the  free-stream 
turbulence  can  not  be  generated  within  the  computational  domain.  To  have  realistic  free-stream  turbulence, 
we  performed  a  large  eddy  simulation  of  isotropic  homogeneous  turbulence  separately.  After  the  turbulent 
flow  evolves  from  its  initial  field  to  a  state  with  specified  intensity  and  length  scale,  we  use  that  flow  field  as 
the  free  stream  turbulence  and  convect  it  into  the  domain  using  the  boundary  conditions  (36). 

Since  the  isotropic  homogeneous  turbulence  is  pre-computed  in  a  finite  domain,  when  it  is  passed  as 
free-stream  turbulence  through  the  inflow,  after  a  certain  time  the  signal  will  repeat  itself.  To  break  this 
periodicity,  phase  jittering  can  be  applied  to  the  free  stream  turbulence  data.  However,  to  minimize  the 
periodicity  of  the  signal,  and  produce  more  realistic  turbulence  for  the  present  flow  configuration,  we  choose 
the  computational  domain  to  be  a  rectangular  box  which  is  4 D,  4 D  and  ID  in  x ,  y  and  2  direction. 
Accordingly,  the  grid  size  is  128x128x32  so  that  Ax  =  Ay  =  Az,  a  necessary  condition  to  preserve  the 
isotropy  of  the  turbulence.  The  numerical  code  for  the  large  eddy  simulation  is  adapted  from  an  existing 
DNS  code  [12]  by  incorporating  the  same  dynamic  SGS  model.  The  temporal  and  spatial  discretization 
are  fourth  order  Runge-Kutta  and  the  sixth  order  compact  finite  difference  scheme  [13].  Periodic  boundary 
conditions  are  applied  in  all  spatial  directions.  The  initial  condition  is  a  purely  solenoidal  velocity  field  with 
uniform  density  and  temperature  field.  The  initial  three  dimensional  energy  spectrum  is  of  the  form 

E(k)  (X  k4 exp[-2(k / kp)2]  (38) 

where  the  peak  wavenumber  of  the  spectrum  kp  is  taken  as  87r.  A  few  small  wavenumbers  which  represent 
the  long  waves  permissible  in  the  x  and  y  directions  but  not  supported  in  z  are  zeroed  out  initially  to  preserve 
the  turbulence  isotropy.  After  the  turbulence  develops  from  its  initial  state,  the  isotropic  turbulence  field 
with  intensity  ufrm8/uoo  =  0.08  and  integral  scale  L/D  =  0.1  is  taken  as  our  free-stream  turbulence  data. 
The  turbulent  Reynolds  number  at  the  inlet  is  Re l  =  u'rm8L/v  =  80. 

7  Results 

The  computational  grid  in  the  x-y  plane  is  generated  by  using  an  algebraic  multi-surface  method  [14]  which 
gives  desirable  orthogonality  at  the  body  surface  and  at  the  inflow/outflow  boundaries.  The  grid  is  clustered 
towards  the  wall  and  the  leading  edge  but  is  uniform  in  z  since  the  flow  is  assumed  homogeneous  in  the 
spanwise  direction.  Unlike  the  conventional  C-type  grid  which  is  largely  stretched  toward  the  far  upstream, 
the  grid  resolution  at  the  inflow  should  be  determined  by  the  characteristic  scale  of  the  free-stream  turbulence 
and  in  the  present  case  should  be  as  isotropic  as  possible.  The  grid  used  in  the  present  simulation  is  of  size 
383x192x64  in  the  streamwise,  normal  and  spanwise  directions.  The  grid  line  distribution  in  x-y  plan  is 
shown  in  Fig. 2.  The  grid  resolution  at  the  inflow  boundary  is  approximately  the  same  as  in  the  LES  of 
isotropic  homogeneous  turbulence.  To  represent  the  experimental  conditions  accurately  with  the  current 
flow  configuration,  the  simulations  are  performed  in  a  progressive  way.  We  first  generate  the  incompressible 
potential  flow  around  the  body  with  blockage  effect  taken  into  account,  then  use  it  as  the  initial  condition  to 
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Figure  2:  Computational  mesh  for  the  upper  half  of  the  domain  in  x  —  y  plane.  For  clarity,  only  every  fourth 
grid  line  of  a  grid  384x192  grid  is  shown  here. 


Figure  3:  Leading  edge  velocity  gradients.  Solid  line  with  plus  sign  is  the  present  simulation,  and  the  circles 
are  taken  from  Van  Fossen  et  al  (1994). 


compute  the  compressible  potential  flow.  After  the  compressible  potential  solution  is  obtained  ,  it  serves  as 
the  initial  condition  and  part  of  the  boundary  conditions  for  the  two  dimensional  laminar  viscous  calculation. 
For  the  three  dimensional  case,  we  take  the  two  dimensional  viscous  solution  as  the  initial  condition  and 
performed  the  simulations  with  the  inlet  disturbances  as  i)  laminar  organized  disturbances  and  ii)  free-stream 
turbulence.  The  purpose  of  studying  the  steady  organized  disturbance  is  to  gain  insight  into  the  basic 
interaction  process  at  the  leading  edge  and  guide  the  choice  of  appropriate  parameters  for  the  free-stream 
turbulence  case.  In  Van  Fossen’s  experiments,  four  different  models  with  the  same  leading  edge  diameter 
but  different  major  to  minor  axes  ratios  from  1:1,  1.5:1,  2.25:1  to  3:1  were  used  to  produce  different  leading 
edge  velocity  gradient.  The  same  model  parameters  are  used  here  for  the  two  dimensional  calculations.  Fig. 
3  shows  the  comparison  of  leading  edge  velocity  gradients  from  the  compressible  potential  solution  M  =  0.2 
with  the  previous  results  obtained  by  a  panel  method[5].  Note  the  leading  edge  velocity  gradients  are  affected 
significantly  by  the  blockage  effect,  and  the  excellent  agreement  shows  that  this  important  feature  of  the 
experiments  has  been  captured  accurately. 

The  two  dimensional  viscous  calculation  is  validated  against  analytical  self-similarity  solutions  at  the 
leading  edge  for  heat  conducting,  compressible  flow  by  Reshotko  and  Beckwith [15].  Fig.  4  and  Fig.  5  show 
the  comparisons  between  numerical  and  similarity  solutions  for  the  velocity  and  enthalpy  profiles  inside  the 
boundary  in  the  leading  edge  stagnation  point  region.  Typically,  the  flow  impinging  on  gas  turbine  blade 
is  in  the  high  subsonic  or  transonic  regime,  so  the  profiles  shown  here  are  for  a  high  Mach  number  case 
M  =  0.8.  The  low  Mach  number  case  M  =  0.2,  which  corresponds  to  the  experimental  conditions,  are  also 
in  excellent  agreement  with  the  similarity  solutions.  In  these  cases  the  ratio  between  wall  temperature  and 
free  stream  total  temperature  are  2.0  and  Pr  is  taken  as  1.0.  As  required  to  possess  similarity  solutions,  the 
viscosity  is  assumed  to  be  a  linear  function  of  temperature. 

Before  introducing  the  free  stream  turbulence  at  the  inflow  for  large  eddy  simulation,  the  inflow  boundary 
conditions  based  Riemann  variants  are  tested.  Fig.  6  show  the  contours  of  the  disturbance  velocities  u, 
v  and  the  vorticity  uj  when  a  circular  Taylor  vortex  is  entering  domain.  The  velocity  disturbance  level 
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Figure  4:  Comparison  of  boundary  layer  velocity  profile  at  the  leading  edge.  Re  =  105,  M  =  0.8,  Tw/To  = 
2.0,  Pr  =  1.0. 


Figure  5:  Comparison  of  boundary  layer  enthalpy  profile  at  the  leading  edge. Re  =  105,  M  =  0.8,  Tw/Tq=2.0, 
Pr  =  1.0. 
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Vmaz/uo c  =  15%  and  radius  of  the  vortex  is  Rrayor/D  =  1.0.  It  can  be  seen  that  vortex  structures  are  kept 
fairly  well  even  for  this  relatively  large  amplitude  disturbance  case. 

The  typical  flow  parameters  for  viscous  laminar  calculations  are  Reo  —  20,000,  =  0.2,  Tw/Tq  — 

1.075.  where  Tw  and  T0  are  the  wall  and  total  free  stream  temperature. The  non-dimensional  leading-edge 
heat  transfer  coefficient,  or  Frossling  number,  is  defined  as 


Fr  = 


Nu 


Re 


1/2 

D 


Tn 

A  TRel2 


(39) 


where  the  Tn  is  the  nondimensional  normal  derivative  of  the  temperature  on  the  body  surface.  As  in  [7], 
the  characteristic  temperature  difference  AT  for  the  evaluation  of  Fr  is  chosen  to  be  the  difference  between 
the  wall  temperature  and  the  adiabatic  wall  temperature  Tw  —  Taw.  The  adiabatic  temperature  can  be 
approximated  by 

Taw  =  \/Pr  (T0  -  Tedge)  -f  Tedge  (40) 

Here  Tedge  is  the  temperature  at  the  edge  of  the  boundary  layer.  Using  the  isentropic  relation  and  assuming 
constant  pressure  across  the  boundary  layer,  the  Tedge  is  obtained  by 


Tedge=Tt( (41) 
To 

where  P0  is  the  free  stream  total  pressure.  For  the  two  dimensional  computation  in  the  absence  of  inflow 
disturbances,  Fr  is  obtained  as  Fr  =  0.748.  The  corresponding  experimental  value  for  the  3:1  model  gives 
Fr  =  0.750  at  a  slightly  higher  Reynolds  number. 

For  the  organized  disturbance  case,  a  spanwise  variation  of  the  form  u  —  uo(x,y)  +  Ausm(2nz/\)  is 
imposed  at  the  inflow  for  the  streamwise  velocity  but  the  total  temperature  at  the  far  upstream  is  kept 
constant.  The  wavelength  is  chosen  to  be  in  the  range  of  0.05D  to  1.0D  and  the  amplitude  is  A u  —  0.05ito- 
Fig.7  shows  the  comparison  of  the  heat  flux  distribution  on  the  body  surface  with  and  without  inflow 
disturbances.  The  value  for  the  disturbed  case  is  averaged  over  the  spanwise  direction.  The  heat  transfer 
rate  is  increased  by  about  30%  with  5%  intensity  of  the  incoming  disturbances  at  A  =  0.1  D.  The  distribution 
of  the  heat  transfer  coefficient  on  the  surface  remains  almost  the  same  as  for  the  undisturbed  case,  which  is 
consistent  with  previous  numerical  and  experimental  studies.  The  dependence  of  heat  transfer  enhancement 
on  the  spanwise  wavelength  A  is  shown  in  Fig.  8.  The  enhancement  decreases  when  the  wavelength  increases. 
But  for  very  short  wavelength,  the  disturbances  are  dissipated  by  viscosity  and  the  flow  reverts  back  to  the  two 
dimensional  case.  So  an  optimum  wave  length  Aopt  which  gives  maximum  heat  transfer  enhancement  is  found 
around  0.1D  for  the  present  flow  condition.  Note  that  even  though  the  stagnation  point  has  the  maximum 
heat  transfer  rate,  the  maximum  increase  m  the  local  Fr  occurs  downstream  of  the  stagnation  point.  When 
large  heat  transfer  enhancement  occuis,  strongly  3D  mushroom-like  flow  structure  are  observed  due  to  the 
vortex  stretching  and  the  large  reversed  streamwise  velocity  at  the  leading  edge.  The  typical  temperature 
contour  projected  onto  the  stagnation  plane  corresponding  to  the  high  heat  transfer  rate  is  shown  in  Fig.  9. 


Before  we  perform  the  LES  of  isotropic  homogeneous  turbulence  in  a  rectangular  box,  the  implementation 
of  dynamic  SGS  stress  model  with  compact  difference  scheme  and  Simpons  3-point  filter  are  first  validated 
through  comparison  with  the  DNS  [16]  and  LES  of  the  Comte-Bellot  and  Corrsin’s  experiments  [17]  in  a 
cubic  box.  Excellent  results  have  been  obtained  which  validate  the  numerical  scheme.  These  are  not  shown 
here  for  brevity.  When  applying  to  the  current  simulation  in  a  rectangular  box,  initially  the  turbulence 
Reynolds  number  is  set  as  ReL  =  100  with  specified  spectrum  and  random  phase  in  spectral  space.  After 
small  scale  eddies  develop  and  the  turbulence  become  realistic,  we  stop  the  simulation  at  ReL  —  80  and  save 
the  data  for  later  use  as  free  stream  turbulence.  Turbulence  integral  scale  L  is  determined  by  a  least-squares 
curve-fitting  to  the  autocorrelation  data  using  an  exponential  function  [5] 

R(r)  =  e~r/L  (42) 

Data  between  0.33  <  R(r)  <1.0  were  used  for  curve  fit.  The  exponential  function  does  not  have  the  correct 
limiting  behavior  for  very  small  values  of  r,  but  the  fit  is  satisfactory  over  the  main  domain  and  the  fitted 
R(r)  can  be  integrated  from  0  to  oo  to  give  the  turbulence  integral  length  scale.  Fig.  10  shows  the  initial 
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Figure  6:  Passing  a  Taylor  vortex  through  the  inflow,  upper  row  t=1.18,  lower  row  t=2.35.  a)  strearawise 
velocity,  b)  transverse  velocity,  c)  spanwise  vorticity. 


Figure  7:  Heat  transfer  rate  on  the  surface  2D  laminar  flow  vs.  spanwise  averaged  3D  case. 


Figure  8:  Increase  of  heat  transfer  vs  spanwise  wavelength  at  different  streamwise  locations  on  the  body 
surface. 
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Figure  9:  Temperature  contours  projected  on  the  stagnation  x-z  plane. 

and  final  three  dimensional  energy  spectrum.  Fig.  11  shows  the  velocity  skewness  reaches  the  typical  value 
~0.4  -0.5  for  the  realistic  turbulence  at  Rgl  about  80.  The  time  development  of  turbulence  kinetic  energy 
is  shown  in  Fig.  12.  Notice  that  wrl  is  slightly  larger  than  ur 2  and  vr2  in  the  time  series,  which  may  be  caused 
by  the  different  size  of  the  computational  domain  in  the  z  direction  compared  to  the  x  and  y  directions. 
Since  this  lack  of  isotropy  is  small,  it  is  not  expected  to  have  any  major  effect.  As  we  started  with  the 
uniform  density  field,  Fig.  13  shows  the  time  development  of  the  RMS  value  of  density  fluctuation.  The 
time  development  of  the  dynamic  model  coefficient  C  is  shown  in  Fig.  14. 

Fig.  15  shows  the  turbulence  intensity  along  the  the  stagnation  stream  line.  The  root-mean-square  values 
are  obtained  by  averaging  u\  v '  and  w*  in  time  and  in  the  spanwise  direction.  The  turbulence  is  largely 
decaying  until  it  reaches  a  distance  of  about  D  from  the  leading  edge  where  the  behavior  of  uf ,  v '  and  w ' 
changes.  Notice  that  close  to  the  body  u '  and  w '  are  amplified  while  v'  continues  to  decay.  In  Van  Fossen’s 
experiment,  a  power  law  curve  of  the  form  Tu  ~  xm  was  used  to  fit  the  decay  of  free-stream  turbulence  in 
the  absence  of  the  model  versus  the  distance  downstream  of  the  turbulence-generating  grid.  The  power-law- 
fitted  curve  is  also  plotted  in  Fig.  15  with  the  same  exponent  m  =  —0.83  as  used  in  the  experiments.  The 
fairly  good  agreement  indicates  that  the  initial  decay  rate  of  the  free-stream  turbulence  has  been  captured 
correctly  by  the  present  simulation. 

Fig.  16  uses  a  logarithmic  scale  on  a  slightly  shifted  x-axis  to  show  the  same  data  as  in  Fig.  15  with  the 
leading  edge  region  amplified  for  clarity.  The  transformation  of  the  x-axis  used  here  is  log( 4.02  -  x).  Note 
x  ~  4  corresponds  to  the  location  of  the  leading  edge.  It  reveals  that  the  free-stream  turbulence  will  go 
through  three  different  processes  as  it  impinges  on  the  leading  edge.  First  is  the  free  turbulence  decay  process, 
corresponding  to  the  distance  from  the  inflow  to  approximately  x  =  3,  when  the  presence  of  the  body  does 
not  affect  the  turbulence  significantly.  The  second  process  is  the  inviscid  distortion  process  approximately 
corresponding  to  the  distance  from  x  =  3  to  x  —  3.95,  where  the  free-stream  turbulence  Uiidergoes  large 
mean  flow  strain  caused  by  the  diverging  mean  streamlines.  A  direct,  quantitative  comparison  between  the 
present  results  and  the  compressible  rapid  distortion  theory  (RDT)[18]  is  not  easily  obtained  due  to  the 
viscous  effect  and  the  relatively  small  scale  of  turbulence,  but  the  increase  of  v!rm8  and  w'rm8,  and  decrease  of 
v'rms  are  qualitatively  in  agreement  with  the  temporal  prediction  of  RDT  under  plane  strain[19].  The  third 
is  the  viscous  interaction,  which  occurs  within  the  leading  edge  boundary  layer  thickness.  Also  notable  is 
that  the  spanwise  velocity  iu'm#  continues  to  increase  significantly  while  u'rm8  and  vlrm9  start  decreasing  due 
to  the  presence  of  the  wall.  This  is  due  to  wall-blocking  or  splatting  effect [20]  [21]. 

Fig.  17  shows  the  vorticity  components  along  the  stagnation  streamline.  Clearly  the  transverse  vorticity 
uy  increases  dramatically  as  the  leading  edge  is  approached.  After  being  wrapped  along  the  surface  as  the 
flow  impinges  on  the  leading  edge,  it  forms  large  streamwise  vorticity  and  produces  the  ’mushroom-like’  flow 
structures  downstream  of  stagnation  point  as  clearly  shown  in  Fig.  21,  which  are  directly  responsible  for 
the  heat  transfer  enhancement.  Fig.  18  shows  the  time  history  of  wrm8  with  the  average  taken  only  in  the 
spanwise  direction.  At  the  inflow  region,  the  presence  of  small  scale  turbulence  is  evident  in  the  fluctuations 
of  the  signal.  By  the  time  the  flow  reaches  the  leading  edge  region,  most  of  the  small  scales  of  turbulence 
are  damped  by  the  viscosity.  However,  under  the  strong  strain  imposed  by  the  mean  flow,  the  free-steam 
turbulence  is  re-organized  into  persistent  flow  structures,  which  significantly  influence  the  boundary  layer  at 
the  leading  edge. 

Fig.  19  shows  the  spanwise  correlation  length  Lc  for  u'  along  the  stagnation  streamline.  Here  Lc  is  defined 
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Figure  10:  Initial  and  final  three  dimensional  energy  spectrum  for  isotropic  homogeneous  turbulence  in  a 
rectangular  box. 


as  one  half  of  the  length  at  which  the  spanwise  autocorrelation  curve  for  uf  first  crosses  zero.  Lc  can  be 
regarded  as  a  length  scale  characterizing  the  size  of  the  energetic  eddies  in  the  spanwise  direction.  A  very 
rapid  decrease  of  Lc  occurs  in  the  leading  edge  region,  indicating  the  strong  vortex  stretching  effect.  Lc 
reaches  its  minimum  value  at  the  outer  edge  of  the  boundary  layer  and  then  increases  inside  the  boundary 
layer  where  viscosity  is  predominant.  The  minimum  value  of  Lc  is  only  about  one  third  of  its  value  at 
inlet.  The  ratio  of  the  minimum  Lc  to  the  local  boundary  layer  thickness  6  is  about  2-3.  These  observations 
suggest  that  the  characteristic  turbulence  scale  at  the  leading  edge  is  largely  determined  by  the  local  strain 
rate  rather  than  by  the  free-stream  turbulence  scale  at  the  inlet.  Computationally,  this  will  determine  the 
required  spanwise  grid  resolution  in  the  leading  edge  region,  because  even  with  large  eddy  simulation,  the 
grid  has  to  be  fine  enough  to  resolve  these  “large  eddies”,  especially  in  high  Reynolds  number  case  where 
the  ratio  of  the  turbulence  length  scales  in  the  free-stream  and  at  the  leading  edge  is  expected  to  increase. 

The  comparison  of  the  Frossling  number  on  the  surface  with  and  without  free-stream  turbulence  is 
shown  in  Fig.  20.  For  the  turbulent  case,  Fr  is  averaged  over  time  and  along  the  spanwise  direction.  The 
increase  in  Fr  is  about  11  %,  and  its  distribution  is  very  much  the  same  as  in  laminar  case.  For  the  same 
flow  configuration,  the  experiments  [5]  were  carried  out  at  a  higher  Reynolds  number.  If  we  extrapolate 
Frossling  number  increase  over  the  laminar  flow  case  using  a  correlation  developed  from  high  Reynolds 
number  experimental  data[5]  to  the  conditions  of  the  current  simulations,  a  value  of  10%  is  obtained.  This 
extrapolation  includes  the  effect  associated  with  the  decay  of  free-stream  turbulence  from  the  inlet  to  a 
location  just  outside  the  boundary  layer*,  as  also  done  in  reducing  the  experimental  data. 

The  current  computational  value  is  in  good  agreement  with  the  experimental  correlation.  The  heat- 
transfer  increase  is,  however,  modest.  This  is  because  of  the  modest  Reynolds  number  of  the  present 
LES.  As  discussed  earlier  the  spanwise  length  scale  decreases  sharply  within  the  boundary  layer,  and  a 
new  computation  with  refined  resolution  is  currently  underway.  Fig.  22  shows  an  instantaneous  heat  flux 
contour  on  the  blade  surface.  Clearly,  the  disturbances  are  largely  stretched  in  the  streamwise  direction,  and 
the  strong  streamwise  vortices  produce  streak-like  regions  with  high  heat  transfer  rate. 


8  Summary 

The  leading  edge  heat  transfer  enhancement  under  incoming  organized  disturbance  and  free-stream  turbu¬ 
lence  are  investigated  numerically.  For  the  laminar  organized  disturbances,  the  three  dimensional  compress¬ 
ible  Navier-Stokes  equations  are  solved  using  high-order  finite  differencing  and  an  implicit  time  marching 
scheme.  Significant  heat  transfer  increase  is  obtained  for  small  inflow  disturbance  amplitude.  The  enhance¬ 
ment  decreases  with  increased  length  scale,  and  an  optimum  length  scale  is  found  which  gives  maximum  heat 

*In  a  conference  paper  AIAA-2001-1016,  we  had  missed  to  account  for  this  change,  and  incorrectly  deduced  that  the  heat- 
transfer  augmentation  observed  in  LES  was  lower  than  the  experimental  correlation. 
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Figure  11:  Time  development  of  the  velocity  gradient  skewness.r  =  L/urm5,  where  L  is  the  turbulence 
integral  length  scale. 
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Figure  12:  Time  development  of  the  turbulence  kinetic  energy. 


Figure  13:  Time  development  of  the  density  fluctuation. 


Figure  14:  Time  development  of  the  dynamic  SGS  model  coefficient  C. 


Figure  15:  Turbulence  intensity  along  the  stagnation  streamline.  The  solid  line  is  a  power  law  fit  using  the 
same  the  exponent  (-0.83)  as  in  Van  Fc^sen’s  experiments. 


Figure  16:  Turbulence  intensity  along  the  stagnation  streamline,  transformed  x  axis. 
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Figure  17:  Turbulence  vorticity  components  along  the  stagnation  streamline,  transformed  x  axis 


Figure  18:  Time  history  of  <  wrl  >1/2  at  locations  close  to  inlet  and  leading  edge  along  stagnation  streamline. 
<>  denotes  the  average  in  spanwise  direction. 


Figure  19:  Spanwise  correlation  length  for  u '  along  the  stagnation  streamline. 
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Figure  20:  Frossling  number  distributions  on  the  leading  edge  region. 


Figure  21:  Cross  section  temperature  contour  at  various  streamwise  locations. 


Figure  22:  Heat  flux  contour  on  the  leading  edge  with  free-stream  turbulence. 
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transfer  increase.  For  the  free-stream  turbulence  case,  a  large  eddy  simulation  using  the  same  algorithm  with 
dynamic  SGS  model  was  performed.  The  free-stream  turbulence  used  in  the  LES  was  generated  by  a  sepa¬ 
rate  large  eddy  simulation  of  homogeneous  isotropic  turbulence.  The  interaction  of  free-stream  turbulence 
impinging  on  the  leading  edge  was  found  to  involve  three  distinct  processes:  free  turbulence  decay,  inviscid 
distortion  and  viscous  interaction.  Strong  streamwise  vortices  are  formed  on  the  wall  through  this  interac¬ 
tion  and  are  largely  responsible  for  the  heat  transfer  enhancement.  The  wall  heat  flux  distribution  under 
free-stream  turbulence  is  found  to  keep  almost  the  same  shape  as  in  the  laminar  case.  Through  the  present 
simulation,  it  is  observed  that  the  typical  energetic  mushroom-like  vortices  which  are  formed  at  the  leading 
edge  are  the  direct  cause  of  turbulent  heat  transfer  enhancement.  The  size  of  these  structures  is  typically 
2-3  times  the  leading  edge  boundary  layer  thickness.  To  characterize  the  scale  of  these  structures  at  different 
Reynolds  number,  as  indicated  in  [8],  the  local  boundary  layer  thickness  would  be  a  better  choice  than  a 
fixed  reference  length  such  as  D .  This  poses  a  severe  challenge  to  a  computational  approach  such  as  large 
eddy  simulation,  because  it  requires  that  these  “large”  energetic  turbulence  structures  at  the  leading  edge 
must  be  fully  resolved  to  have  an  accurate  prediction  of  turbulent  heat  transfer  under  free-stream  turbulence. 
For  the  modest  Reynolds  number  of  the  current  simulations  the  predictions  were  found  to  agree  well  with 
experimental  data.  However,  to  capture  the  free-stream  turbulence  effects  in  higher  Reynolds  number  flows 
it  will  be  necessary  to  use  a  significantly  refined  grid  within  the  boundary  layer.  Our  results  suggest  that 
what  is  needed  is,  in  effect,  a  hybrid  simulation  strategy  where  the  turbulence  outside  and  away  from  the 
blade  surface  is  captured  using  LES  techniques  while  a  finer  DNS-like  grid  is  embedded  within  the  near- wall 
region  to  resolve  the  smaller  eddies  responsible  for  near-wall  effects.  Such  a  methodology  is  being  developed 
in  an  extension  of  the  work  supported  under  this  grant. 

One  of  the  primary  objectives  of  current  research  is  to  develop  large  eddy  simulations  as  a  tool  for  heat 
transfer  prediction  over  a  turbine  blade  immersed  in  a  high-subsonic/ transonic  hot  stream  containing  free 
stream  turbulence.  So  the  next  logical  step  wall  be  to  extend  the  current  simulation  to  a  higher  Mach 
number  case,  which  is  not  easily  investigated  experimentally.  The  three  dimensional  moving  wake  caused 
by  the  upstream  blade  rows  is  another  essential  element  of  the  flow  impinging  on  a  downstream  blade  row. 
Instead  of  isotropic  homogeneous  turbulence,  wake  turbulence  can  also  be  introduced  as  the  free-stream 
disturbance.  This  combined  study  would  provide  a  unique  assessment  of  the  important  blade  heat  transfer 
mechanism  relevant  to  the  gas  turbine  environment. 
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